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1. INTRODUCTION 

Maximum power point tracking (MPPT) is of the paramount importance in renewable energy 
regimes for not only to maximize the system’s power efficiency, but also to reduce the return period of the 
installation expenses cost and good power quality and reliable. In this research, the targeted algorithm is hill 
climbing search (HCS) under the direct power control (DPC), to increase the output power direct through the 
duty cycle, and to reach maximum peak point on the DC link. HCS method is considered an example of a 
perturb and observe (P&O) technique [1-4]. There are many features for HCS algorithm; it’s considered the 
simplest method, not require any prior knowledge about wind energy conversion system (WECS) or 
aerodynamics characteristics. It can be applied to any WECS, suitable for small scale wind turbine (WT) and 
this type of algorithm can bring the operating point toward power coefficient (Cp) by increasing or 
decreasing the perturbing [5-9]. 

There are two main types according to the maximization of captured power. The first type is DPC 
concern in the power output direct in the WECS. Whereas, the second type is indirect power control (IDC), 
increased the mechanical power in the WECS as shown in Figure 1. 

The input of HCS is considered electric power since it can be simply measured using the converter 
[10-15]. However, it is in reality the turbine power that must be applied for the control strategy to reach to the 
peak point. Since no data for the wind turbine WT is needed, it is ensured that the WT will reach to its real 
maximum power point, even through variations of external blade characteristics or factors [16]. 
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In Figure 2, HCS algorithm concept, perturbing the control variable until reach maximum peak 
point. In spite of these special properties that lead HCS algorithm to be the best option for MPPT control in 
any WECS. However, in fact, it is only appropriate in the slow changing wind speed conditions. Therefore, 
HCS traditional algorithm has two major problems; trade-off between speed and efficiency, and the second is 
the perturbation step size [8, 17-19]. 
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Figure 1. Classification of MPPT according the maximization 
of total captured power 
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Figure 2. The main principle of HCS 
algorithm, Perturbing the control variable 
until reach maximum peak point 


The main objective of this research is to increase the power efficiency for the WECS. The power 
efficiency for the WECS can be increased using HCS algorithms to decrease the number of iterations to reach 
the decision for the MPPT, and tuning parameters for the DC link [20]. Palestine is suffering from 
environment pollution problems, especially Gaza strip region after the last three wars. In addition, the siege 
imposed on Gaza strip, and the continued interruption of electric power as well as fuel make the need to an 
alternative source of energy instead of traditional sources. The Gaza Strip has an increasing demand for 
electrical power with an increasing shortage of power supplies [21-23]. Figure 3 shows Palestine map which 
there is a deficit in electricity due to siege and prevailing war-torn conditions. Moreover, electricity 
generation is not feasible to satisfy the demand. 

Many researchers are actively proposing algorithms to extract maximum power point out of WT. 
However, this solution cannot be adaptable for Palestine due to the different wind condition compared with 
other places [9]. There are many techniques under MPPT. However, among these, three mode HSC algorithm 
is more feasible to be implemented in Palestine, this is due to reduction of step fluctuations and suitable for 
micro grid system, but the three mode HCS algorithm cannot reach maximum powerpoint due to 
convergence speed delay, wrong tracking direction, trade-off between speed for detection the maximum 
power point and accuracy of tracking technique [26-28]. 

Figure 4 shows the three mode HCS algorithm [17]. Many systems have a master controller that 
knows in which mode the controller is operating, relaying on wind speed, differences in wind speed, in this 
approach, the controller can response differently on smaller or larger wind speed changes or it can keep the 
rotor speed steady as long as the wind speed changes does not exceed any dead band limit. The probability 
for adaptive HCS algorithm arises, which are first operating in a learning mode to determine all the important 
parameters in function of a certain wind pattern [29]. 

Figure 5 illustrates the mechanism for the three mode HCS algorithms to reach the MPPT. 
Basically, it’s relay on the comparison between the pervious state and the next state. Thus, it spends long 
time to reach the maximum peak point. Moreover, the algorithm is very complicated due to the iteration 
process to reach MPPT [8, 17]. 
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Figure 3. Palestine map [24, 25] 
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Figure 4. Principle of the HCS three mode 
algorithm [17] 





Figure 5. Flowchart of three mode HCS algorithm MPPT [17] 
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2. WIND SPEED AND MPPT RELATION 

The power output from a WT relies on the rotor swept area, density of air and velocity of the wind. 
The common power equation for a WT is given by (14) based on case of constant acceleration [30]. The 
velocity v is equal to the work done W, and the kinetic energy of an object having mass m, the distance s 
under a force F, as the following. 


E =W = Fs (1) 
Applying Newton’s Law: 

F = ma (2) 
Hence, 

E = mas (3) 
Using equation (3) of motion: 

v? = u? + 2as (4) 
We can get: 


_ (v? —u*) 


2s 


(5) 
The object initial velocity is 0, i.e. u = 0, then: 

a = — (6) 
Substituting it in (3), the kinetic energy of a mass is: 

E = 0.5mv? (7) 


The wind power is given by the rate of change of energy as the following equation: 


dE 1 d 
P = =? = (8) 
dt 2 dt 


As mass flow rate is given by: 


dm ax 

a Pa (9) 
The rate of change of distance is given by: 

ax 

a) (10) 
Then 

dm 

Fos pAv (11) 


E=Kinetic energy (J), p = Density (kg/m°), m=Mass (kg), A=Swept Area (m^), v=Wind Speed (m/s), 
Cp=Power coefficient, P=Power (W), dt/dm=Mass flow rate (kg/s), r=Radius (m), x=distance (m), t=time (s), 
dt/dE=Energy flow rate (J/s). 


Therefore, from (8), the power can be represented as the following equation: 


P =0.5pAv? (12) 
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The captured power from the wind can be expressed as the following equation: 
Pavail = 0.5pAv?Cp (13) 


The swept area of the wind turbine can be computed from the length of the wind turbine blades 
using the equation for the swept area of acircle A = mr 


P= 0.5cpa g PTR? Vy (14) 


Where, cp = power coefficient, Vw = Wind velocity, pa = density of air, and and R = rotor radius. The power 
coefficient (cp) for the WT is the maximum power that could be captured of the total available in the wind, 
cp is a function of blade pitch angle (f) and blade tip speed ratio (4). Theoretically, a wind turbine can extract 
maximum 59.3 % from the total power of wind (Betz’s limit)[9]. Betz limit identified as a WT can only 
convert 59.3 % maximum of the total available energy in the wind into kinetic energy. Practically a WT can 
extracts around 40 % of the power in the wind in the ideal WT power [31]. 

For small scale WT, p angle between the reference line and chord of the blade on the rotor hub, 
normally zero or kept constant. Consequently, for a constant pitch angle, power coefficient is only a function 
of tip speed ratio. The tip speed ratio A is a ratio between the wind speed (Vw) and speed of the blade tip 
(Vtip) as shown in (15) [31]. 


— ip. “Or 
Nee ae (15) 


where, œr is the turbine angular speed. The wind turbine dynamic equation is given by (16) [31]. 


dwr = Tm—-TL-Fw,y 
a (16) 





where, F is the viscous friction coefficient, J is the system inertia, TL is the torque due to load, Tm is the 
torque given by the turbine. The power coefficient Cp(A, p) can be determined for WT by (17). 
= Ce 


Cy (A, f) = ag — c3ß — c4)e ĉi cel (17) 


where, Ai is a variable which is a function of J, p is defined as in (18), 


1 1 0.035 
Ai A+0.08B ßB3=+1 (18) 

Substituting (16) through to (17), the coefficients from c1 to c6 are: c3 = 0.4,c1 = 0.5176, c4 = 5, c2 
= 116,c6 = 0.0068 and c5 = 21 [31]. Figure 6 represents the nonlinear power coefficient cpa, p, described by 
(17), according to the tip speed ratio À for various values of B pitch angle [31]. 

Mean wind speed play important role in the wind power generation process due to the cubic 
proportional relation between power and mean wind speed. Average wind speed reflects to the rotor speed. 
Therefore, total output power is proportional with rotor speed as shown in (14). Figure 7 presents that rotor 
speed versus output power for each wind speed, there is a maximum output amount of power that the WT 
could extract if WT operated at the rated rotor speed (@opt) or to the optimum rotor speed value Wo, can 
reached at the optimum (A,,;). in order to collect maximum possible power from the WT. WT must be turned 
on at optimal TSR 4,,;. This case of operation is possible by controlling the rotational speed of the WT that it 
always turns on at the optimum or rated speed. 


Prax = S (19) 
Where, Kopt 18 an optimum wind written by (20). 
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3. PROPOSED HCS ALGORITHM TO REACH MPPT 

In this work, three mode algorithms had been combined to be simpler. The adaptive algorithm 
contains three mode, in it’s considered very complicated due to iteration for same equation of w*[8, 33]. 
Figure 8 shows the novel two mode proposed algorithm. The main aim from the new algorithm is to increase 
power efficiency, increasing the convergence speed and avoid the trade-off between speed for detection the 
maximum power point and accuracy of tracking technique. The decision of the algorithm through duty cycle 
can reach to the curve to maximum peak point due to fewer calculations in the mode stages. The new 
algorithm is considered simpler than the three modes and faster. 
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Figure 8. Flowchart of the proposed novel HCS algorithm. 
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4. MODELING THE WECS USING PMSG WIND GENERATORE 

There are five main blocks in the WECS model, which are WT, permanent magnet synchronous 
generator (PMSG), three phase diode bridge rectifier, controller and DC to DC boost converter. Six wind 
speed profile effect to the WT as input. Scopes in PMSG measured the stator current, rotor speed, 
electromagnetic torque and stator voltage. The three-phase output from the PMSG has been rectified using 
the three-phase diode bridge rectifier. 

Figure 9 illustrates the WECS model with three mode HCS algorithm and the proposed HCS 
algorithm. In this model it can be applied the HCS algorithms to check the efficiency performancefor each 
algorithm based on the output signal behavior. The scope measurement has been checked the DC voltage via 
DC link. To enhance the MPPT technique it should be tuning for the model circuit especially DC to DC boost 
converter. The main function of the controller is to make switching through the duty cycle based on the 
decision of the algorithm. The three-mode algorithm (Adaptive) makes iterations to reach the MPPT.whereas, 
the proposed algorithm avoids the iteration delay time to reach the MPPT due to the compined mode as 
shown in Figure 8. 
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Figure 9. WECS-PMSG model 


4.1. WECS PMSG with three mode algorithm (adaptive) 

In this section, the adaptive or the three-mode algorithm has been applied to the WECS PMSG 
model using MATLAB SIMULINK. Six wind speed profile the input of the WT. The MATLAB scope 
measured rotor speed, stator voltage and current, electromagnetic torque and power for PMSG. To enhance 
the results for the PMSG it should be tuning the parameters to reach to the targeted objective. 

In Figure 10 DC voltage had been measured before applying the three-mode algorithm to describe 
the signal behavior before and after applying the three-mode algorithm, after that the effect of the three-mode 
algorithm and the proposed algorithm 1s compared. 

Figure 11 shows the DC voltage after applying three mode algorithms. There is a delay to reach the 
maximum power point because of an iteration processes to reach the MPPT. The three-mode algorithm had 
tradeoff between efficiency and convergence speed to reach max peak point due to complexity to reach the 
decision for the peak point. 
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Figure 10. DC voltage before applying the three-mode algorithm 
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Figure 11. DC voltage after applying HCS three mode algorithm 


4.2. Appling the proposed novel algorithm 

The proposed algorithm has been applied to the WECS. Six wind speed profile the input of the WT 
to check the efficiency performance for the model under fast varying in wind speed. In Figures 12, 13 it can 
be seen the three Phase voltage and current are pure sine wave which amplitude increased dramatically direct 
based on increasing wind speed. Voltage curve and current are in phase, there is no leading or lagging in the 
phase angle which equals zero. Therefore, the power factor equals 1. It means that inductive reactance equals 
to capacitive reactance and that condition called resonance condition. 
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Figure 12 Zoom of generator currents iapc between 4.672s and 4.678s 





Figure 13. Zoom of generator Voltage Vape between 14.99s and 15.07s 


The simulation results are shown in pervious figure, the turbine reach to its maximum peak point at 
a wind speed of 15 m/s is applied. The electric power and rotor speed are sampled every 3 seconds after a 
new rotor reference speed is set. The frequency of the PMSG generator is equal to 50 Hz. 
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The new algorithm is controlled in D duty cycle to control the voltage to reach maximum peak 
point. Bitch angle is considered constant (B=0) irrespective of wind speed. It’s clear from Figures 15 and 16 
after the rating wind speed the proposed algorithm reached the curve to the maximum point value through the 
duty cycle. Figure 16 shows the performance of the proposed MPPT algorithm under mixed wind profile 
with six variation point. The performance of the developed algorithm presents fast tracking capabilities. 
Moreover, the new technique for HCS can reached to the maximum peak point with minimum calculation 
which making it simple implementation algorithm and competitive due to summarized three mode to two 
mode. The proposed algorithm shows fast tracking capability and enhanced stability under both low and high 
rate of change wind speed conditions and is verified using MATLAB/Simulink. 


Figure 14. Zoom of generator power output between 15.00 s and 15.06s 


Figure 16. Voltage after applying two mode algorithms for WECS-PMSG 


5. CONCLUSION 

This paper proposed a new technique for HCS algorithm to increase the efficiency performance to 
reach to the maximum peak point. The proposed algorithm captured the maximum possible power and 
reached to the maximum power point on the curve at different wind speed without the knowledge of the WT 
aerodynamic characteristics. The performance of the proposed novel algorithm presents fast tracking 
capabilities. It can reach to the maximum peak point with minimum calculation which making it simple 
implementation algorithm and competitive due to summarized three mode to two mode. The experimental 
results confirm that the proposed technique for HCS is remarkably faster by 20% of the total time required 
comparing to the mode algorithm due to simplicity. 
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